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ABSTRACT 

We compare the gravitational potential profiles of the elliptical galaxies NGC 4486 
(M87) and NGC 1399 (the central galaxy in the Fornax cluster) derived from X-ray 
and optical data. This comparison suggests that the combined contribution of cosmic 
rays, magnetic fields and micro-turbulence to the pressure is '^10% of the gas thermal 
pressure in the cores of NGC 1399 and M87, although the uncertainties in our model 
assumptions (e.g., spherical symmetry) are sufficiently large that the contribution 
could be consistent with zero. In the absence of any other form of non-thermal pressure 
support, these upper bounds translate into upper limits on the magnetic field of ~10- 
20 ^G at a distance of l'-2' from the centers of NGC1399 and M87. We show that 
these results are consistent with the current paradigm of cool cluster cores, based on 
the assumption that AGN activity regulates the thermal state of the gas by injecting 
energy into the intra-clustcr medium. The limit of ^^10-20% on the energy density in 
the form of relativistic protons applies not only to the current state of the gas, but 
essentially to the entire history of the intra-cluster medium, provided that cosmic ray 
protons evolve adiabatically and that their spatial diffusion is suppressed. 
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1 INTRODUCTION 

Both optical and X-ray data are often used to determine the 
distribution of the gravitating mass in galaxies. In analysing 
optical stellar-kinematic data, stars are treated as collision- 
less particles and the mass distribution is obtained either 
from the Jeans equations and some assumption about the 
anisotropy of stellar orbits (Binney & Tremaine 2008; Lokas 
& Mamon 2003), or by finding the potentials in which a dis- 
tribution of orbits reproduces the observed surface bright- 
ness profile and stellar kinematics (e.g., Kronawitter et al. 
2000; Thomas et al. 2007). In X-rays, hydrostatic equilib- 
rium is usually assumed for the gaseous atmosphere and 
the observed density and temperature distributions are then 
used to evaluate the distribution of gravitating mass (Math- 
ews 1978; Forman, Jones & Tucker 1985; Fukazawa et al. 
2006; Humphrey et al. 2006). Ideally both methods should 
yield identical results, and the mismatch between gravitat- 



ing masses derived from optical and X-ray data can be used, 
for example, to measure the contribution of non-thermal 
particles to the gas pressure. Such non-thermal particles 
are directly observed as bubbles of relativistic plasma (e.g., 
Boehringer et al. 1993). They can also be present in the ICM 
due to mixing of thermal and relativistic plasma, or gener- 
ated by shocks propagating through the ICM. Below we use 
the newest Chandra data to compare optical and X-ray data 
in two giant elliptical galaxies. The ~ 0".5 angular resolu- 
tion of Chandra has brought the X-ray data on a par with 
optical data and has made it possible to extend the com- 
parison between the kinematics of stars and hot gas from 
spatial scales of arcseconds up to ~ 10' as is done below for 
two giant elliptical galaxies: M87 and NGC1399. 

The fundamental assumption that underlies this pa- 
per is that the mass profiles derived from optical obser- 
vations are correct, so that deviations between the opti- 
cal and X-ray profiles are due to effects on the gas such as 
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non-thermal pressure. While measuring mass profiles from 
stellar-kinematic data involves only gravity and dynamics, 
uncertainties may remain, mainly because of the degener- 
acy between radial variations in mass-to-light ratio and ra- 
dial variations in the anisotropy parameter of the velocity 
ellipsoid (Binney & Mamon 1982). Breaking this degeneracy 
is possible when the kinematic data provide sufficient infor- 
mation on higher-order moments of the line-of-sight velocity 
distribution (LOSVD) and this is used to constrain the dis- 
tribution of stellar orbits (Gerhard 1993, Merritt 1993). Un- 
certainties resulting from modeling good quality absorption 
line spectra are ~ 10% in circular velocity over the range 
of radii covered (e.g., Kronawitter et al. 2000; Thomas et 
al. 2007, and references therein). Additional reasons to be- 
lieve that the optical mass profiles that we use in this paper 
are reliable are: (i) the profiles from M87 and NGC1399 
are derived from quite diflerent methods — based on stel- 
lar absorption line spectra but including also discrete tracer 
particles (globular clusters) in the former case — yet yield 
consistent results for the level of non-thermal pressure; (ii) 
for both galaxies, the mass distributions have been stud- 
ied in at least two independent studies that yield consis- 
tent results; (iii) most of the concerns about the reliabil- 
ity of optical mass profiles arise from studies with limited 
observational resolution (e.g., determinations of black-hole 
masses) , limited radial coverage of data, or restrictive mod- 
eling techniques (e.g., the assumption of isotropic velocity 
distributions), neither of which is a concern in the studies 
used here. 



2 SELECTION OF M87 AND NGC1399 

The choice of M87 and NGC1399 is motivated by the fol- 
lowing considerations. We want the objects (i) to be nearly 
round (EO-El) elliptical galaxies with existing high quality 
optical data on stellar kinematics, (ii) to have X-ray emis- 
sion dominated by diffuse gas, with temperatures not much 
higher than the velocity dispersion of the galaxy, (iii) to have 
high-quality Chandra data; (iv) to exhibit at most small or 
moderate deviations from spherical symmetry in the X-ray 
images of the hot gas. There are few objects that satisfy all 
these criteria. For example, condition (i) requires that the 
object be nearby, while (ii) requires that we study massive 
early-type galaxies at the center of a group or a cluster with 
a cool core. Choosing a massive galaxy is necessary to avoid 
contamination of the X-ray emission by unresolved point 
sources. For instance in NGC3379 (the dominant elliptical 
galaxy in the poor Leo group), the X-ray emission remaining 
after removing bright low mass X-ray binaries has a large 
(perhaps dominant) contribution from the unresolved stel- 
lar population (CVs, coronally active stars). This conclusion 
follows from the comparison of the unresolved X-ray lumi- 
nosity per unit K-band optical luminosity of NGC3379 with 
that for the gas poor bulge of M32 (Revnivtsev et al., 2007, 
2008; see also David et al. 2006). Application of the hydro- 
static equilibrium equation to such systems would lead to 
the conclusion that the gas has been heated to a tempera- 
ture that cannot be bound to the galaxy, and is out-flowing 
(Pellegrini & Ciotti 2006). 

Both M87 and NGC1399 are at least several times more 
massive than NGC3379. They reside in the gas-rich cool 



cores of nearby clusters (Virgo and Fornax respectively), and 
both objects have excellent optical data and long Chandra 
observations (see Fig |ll2|l . Thus, although we plan to apply 
our technique to additional galaxies in the future, M87 and 
NGC1399 satisfy all criteria enumerated above and are the 
best targets for our pilot project. 

We assume distances of 16 and 19.8 Mpc for M87 and 
NGC 1399 respectively 



3 OBSERVATIONS AND DATA 
PREPARATION 

Our Chandra data analy sis is adopted from th e reduction 
procedure described in (|Vikhlinin et al.l |2005| . V05 here- 
after) . This includes filtering of high background periods and 
application of the latest calibration corrections to the de- 
tected X-ray photons, and determination of the background 
intensity in each observation. 

The quiescent Chandra background is dominated by 
the cosmic X-ray background and charged particle events. 
The latter component can be subt racted exquisitely ac- 
curat ely (with a < 2% scatter, see iHickox fc MarkevitchI 
120061 ) if one removes the secular trends in background in- 
tensity using the hard X-ray band data. The cosmic X-ray 
background component is modeled adequately by using the 
"blank-sky" background dataset^J which include both the 
particle-induced and unresolved sky components. 

The spectral calibration includes removal of the time- 
dependent gain variations of the ACIS CCDs, and also the 
most recent corrections for all positional-dependent varia- 
tions in the ACIS response and effe ctive area. For the discus- 
sion of associated uncertainties, see lVOSl . They are negligible 
for our purposes. 



3.1 M87 data 

For the analysis we used Chandra observations (OBSIDs 
5826, 5827, 5828, 6186, 7210, 7211, and 7212) taken at a 
variety of instrument roll angles from February to November 
2005 using the ACIS-I detector (CCDs 10-13) in Very Faint 
(VF) mode to minimize the background (see Forman et al. 
2007 for details). 

We reprocessed all observations applying the latest 
CTI and time-dependent gain calibrations. We performed 
the usual filtering by grade, excluded bad/hot pixels and 
columns, removed cosmic ray 'afterglows', and applied the 
VF mode filtering. We also reprocessed two Faint mode 
ACIS-S OBSIDs (3717 and 2707) and treated the front (S2) 
and back (S3) illuminated CCDs independently. We com- 
pared the images from the Very Faint and Faint mode ob- 
servations to verify that no artifacts were introduced near 
the bright jet by the use of the large 5x5 pixel event regions 
in VF mode. We examined the data for background fiaring 
and found moderate fiaring in OBSIDs 3717 and 2707 (the 
back illuminated CCDs only) for which we excluded approx- 
imately half the duration. A typical effective exposure time 
is ~ 500 ksec. The background files (see Markevitch 2001 



^ See |http://cxc. harvard. edu/contrlb/m£Lxim| 
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for details) were processed in exactly the same manner as 
the observations. 

During 41 ms of the 3.2 s nominal exposure readout of 
the ACIS CCDs, the chips are exposed to the sky. This re- 
sults in a small contribution of the source flux, 1.3%, being 
uniformly re-distributed along the readout direction. This 
'readout artifact' is most clearly visible when there is a 
bright point source in the FOV which produces a charac- 
teristic readout streak. However, for any emission where a 
bright region contaminates one of lower surface brightness, 
this effect must be taken into account. The readout artifact 
can be accurately subtracted using the technique described 
by Markevitch et al. (2000) . A new event file is generated us- 
ing the original data by randomizing the CHIPY-coordinate 
and all sky coordinates and energies are recalculated as if it 
were a normal observation. The new event file is renormal- 
ized by the ratio of the readout time compared to the inte- 
gration time (41 ms/3.2 s = 0.013) and then treated exactly 
as another component of the background. The blank-field 
background is also renormalized by reducing its integration 
time by 1.3% to account for this additional subtraction. The 
M87 nucleus and the brightest jet knot were piled-up in most 
of the observations and the corresponding readout streaks 
were completely excluded from the analysis. 

As the last step of data preparation, an additional col- 
umn was added to each event list. For a given event in the 
list this column 77 contains the ratio of the effective area for a 
photon with a given energy at a given position A{E, Xd, yd) 
to a predefined function of energy Aq{E): 



rj^A{E,Xd,yd)/Ao{E), 



(1) 



where A{E,Xd,yd) includes mirror and detector efficiencies 
(including non-uniformity of the detector quantum efficiency 
and the time and spatially dependent contamination on 
the optical blocking filter). Our data set largely consists of 
ACIS-I data and a natural choice is to set Ao{E) to the 
ACIS-I on-axis effective area. This makes 77 equal to the vi- 
gnetting of the mirrors, modified by the (energy dependent) 
variations of quantum efficiency across the detector. The 
same procedure was repeated with the background data. Fi- 
nally, for each event list, we make an exposure map that ac- 
counts for all position dependent, but energy independent, 
efficiency variations across the focal plane (e.g., overall chip 
geometry, dead pixels or rows, variation of telescope point- 
ing direction). 



3.2 NGC 1399 data 

For NGC 1399 we used Chandra data of OBSIDs 319, 4172 
and 4174, having a total useful exposure of ~120 ksec. The 
data were cleaned and processed using the same procedures 
as for the M87 data. 
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Figure 3. M87 surface brightness in the 0.6—2 keV band (lower 
panel) and the cmissivity profile in the same band (upper 
panel). The solid line in the upper panel shows the function 
(1 + (r/rc)2)"^'3, with /3 = 0.37 and re = 0.2'. Clearly, the /3- 
niodel provides a reasonable, but not perfect, description of the 
gas cmissivity in the 0.6-2 keV band. The thick green line in 
the lower panel shows the projection of the cmissivity profile for 
comparison with the observed surface brightness. It goes through 
the data points, showing that the deprojection solution is exact. 
The nearly horizontal dotted line is the contribution to the sur- 
face brightness in each annulus from the outer gas layers (outside 
19'), where the gas cmissivity is assumed to decline as a power 
law of radius. 



4.1 Broad band 

First, using the event lists for all 9 observations, we build an 
image in the 0.6-2 keV band, centered at the M87 nucleus, 
40' on a side and with l" pixels. To build this image, the 
image pixels were initialized to zero, then for each event in 
the event list with energy in this range, we added the value 
l/»7 to the corresponding image pixeo. The same proce- 
dure was applied to the background event lists. An exposure 
map was generated combining exposure maps for individual 
event lists. A background-corrected surface brightness pro- 
file Sobs{i) calculated from the '1/r;' image in a set of annuli 
centered at the M87 nucleus is shown in the bottom panel 
of Fig. |3] The surface brightness is normalized to counts per 
arcmin^ per second. 

In our deprojection analysis, we follow the approach 
used in Churazov et al. (2003) (see Appendix |X] for addi- 



4 DEPROJECTION 

We use the observations from M87 to describe our method in 
detail and to explore biases arising from multi-temperature 
gas components. Such components are clearly visible in M87 
(see Fig.[T} and are less prominent in NGC 1399 (see Fig. [2]). 
Hence, M87 provides the most stringent test of our method. 



^ Since r] is essentially the vignetting factor (relative to the on- 
axis position), populating images with l/rj instead of 1 for each 
event is equivalent to applying the vignetting correction during 
image construction instead of accounting for this factor in the 
exposure map. The logic behind this approach is further discussed 
in Section 1X21 
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Figure 1. Central 10' X 10' (10' = 47 kpc) region of the Chandra 0.6-2 keV (left) and DSS optical image (right) of M87. The galaxy 
appearance is very regular in the optical band, while the X-ray image is moderately disturbed. In these and other images throughout 
the paper, north is up and east is to the left. 




Figure 2. Central 10' X 10' (10' = 58 kpc) region of the Chandra 0.6-2 keV (left) and DSS optical image (right) of NGC1399. As for 
M87, the optical galaxy is very regular, while the X-ray image shows only moderate deviations from spherical symmetry and appears 
considerably less disturbed than M87. 



tional discussion) . We assume spherical symmetry, but make 
no specific assumption about the form of the underlying 
gravitational potential. For a given surface brightness pro- 
file in Ua annuli, we choose a set of n^ (n^ ^ Ua) spherical 
shells with the inner radii r(i), i = 1, . . . , n^. The gas emis- 
sivity £ is assumed to be uniform inside each shell, except 
for the outermost shell, where the gas emissivity is assumed 
to decline as a power law of radius: £ = £outr~^'^''"\ where 



Pout is a parameter. In our analysis r{ns) was set to 19' 
and Pout = 1/3- In practice, we are mostly working with the 
data inside a 10' circle and the precise values of Pout and 
the outer radius do not affect the results. 



The expected surface brightness can then be written as 
a projection of emissivities in each of the n^ — 1 shells plus 



a contribution from the outer layers (r > r{ns)): 
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SU) = E ^(^.i)^W + Pout{j)£ouU 



(2) 



where £{i) is the emissivity of a given shell, £out is the emis- 
sivity of the outer layers of the gas at r = 1, and P(i,j) 
and Poutij) are the projection matrix/vector from our set 
of shells into our set of annuli. A simple analytical expres- 
sion for P{i,j), which is a function of the geometry only, is 
given by McLaughlin (1999). 

Deprojection can then be reduced to a simple least 
squares problem - what set of emissivities in our set of 
shells (together with the emissivity normalization £out for 
the outer layers) provides the best description of the ob- 
served surface brightness: 



X 



J2lS{j)~Sots{j)f/aijf 



(3) 



where a{j) is the error associated with the surface bright- 
ness in annulus j. In our analysis, we use modified errors 
'^0)^ = f^statUY +5^ X Soi,s(j''), where Ustat{j) is the statis- 
tical error associated with the Poisson noise in the observed 
image and background and 5 ~ 0.1. This was done to avoid 
the situation where a few annuli with the best statistics com- 
pletely dominate the 'x^ . In practice for the M87 dataset, the 
choice of 5 does not affect the results. As usual the differen- 
tiation of this relation with respect to £{i) and Eout yields 
a system of linear equations P£ — S, which can easily be 
solved. Here P is the square Us x Us matrix and each of the Us 
elements of S vector are the linear combinations of the orig- 
inal observed values of Sobs- In the cases considered below, 
the solution of equation Q is unique, once the parameter 
Pout is fixed. The properties of the inverse matrix P~^ (in 
particular an error enhancement when the ratio of maximal 
to minimal eigenvalues of P is large) can easily be controlled 
by making the spherical shells broader. A practical recipe for 
the choice of shells/annuli radii will be described elsewhere. 
The emissivity of each shell can then be evaluated as 
an explicit linear combination of the observed quantities: 

£{z),£out^ J2 P'\j,i)S{j)^ J2 D{j,i)Sois{j), (4) 



J — IjTls 



J — IjTla 



where D is the final Ua x Us deprojection matrix relating 
the observed surface brightness in rings and the emissivities 
of spherical shells. Since the whole procedure is linear, the 
errors in the observed quantities can be propagated straight- 
forwardly. Shown in Fig. [3] (upper panel) is the set of gas 
emissivities in the 0.6-2 keV band derived from the surface 
brightness profile. The forward convolution of this set of 
emissivities with the projected matrix P{i,j) (green line in 
the lower panel of Fig. |3} exactly reproduces the observed 
surface brightness profile. 



4.2 Deprojected spectra 

Note that with our definition of 77, the projection ma- 
trix P{i,j) does not depend on energy. Indeed all en- 
ergy/position dependent factors are already corrected for 
when we make the I/77 image instead of a usual image con- 
structed from raw counts without any weights. We do lose 
some of the sensitivity with this I/77 procedure, but very 
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Figure 4. Projected (top) and deprojected (bottom) spectra for 
the 7'. 4 — 8'. 8 annulus/shell in M87 (black crosses) together with 
the best fitting single-temperature APEC model (red and blue 
histograms). The best fitting temperatures are 2.57 and 2.4 keV 
for projected and deprojected spectra respectively. Both spec- 
tra are reasonably well described by a single temperature APEC 
model. 



little as long as the variations of l/ry across the region used 
for spectrum extraction are not large. Note that since the 
exposure time is dominated by that from the ACIS-I chips, 
the factor 77 is a smooth function of energy and our l/rj 
technique does not introduce spurious spectral features. 

Since P{i,j) does not depend on energy, the deprojec- 
tion of a surface brightness profile in any energy band would 
look like equation (|4]) with exactly the same D(j,i). Thus 
we can accumulate a set of 'I/??' spectra (corrected for back- 
ground and readout) for each of the Ua annuli, and apply 
equation (|4]) to determine the emissivities of each shell in 
each of the ACIS energy channels. We note here that a simi- 
lar technique of vignetting correction at the level of individ- 
ual events is often used in the analysis of X-ray data (e.g. 
for XMM-Newton) . For M87 Chandra observations we ex- 
plicitly verified that fitting of 'l/jy' and 'regular' spectra of 
the same regions yields fully consistent results. 

The procedure described above differs slightly from the 
often used onion peeling algorithm or pro jet in XSPEC 
(Arnaud 1996). A short note on the different deprojection 
techniques is given in the Appendix |X] The bottom line is 
that all techniques should yield consistent results if under- 
lying assumptions are correct. Our procedure is computa- 
tionally fast and is especially useful when a large number of 
spherical shells are analyzed. 

The resulting spectra were fitted in XSPEC V12 (Ar- 
naud 1996) with a single-temperature optically thin plasma 
emission model as implemented in the APEC code (Smith 
et al. 2001). The low energy photoabsorption was fixed at 
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Figure 5. Radial profiles of deprojected gas parameters in M87. 
The parameters were obtained using a single temperature APEC 
model (with fixed low energy absorption and redshift) fit to the 
deprojected spectra for a set of spherical shells centered on M87. 
The red and blue points correspond to the best-fitting parame- 
ters for the 0.6—9 keV and 2.5-9 keV energy bands respectively. 
Simple analytical approximations for the density and tempera- 
ture are shown with the solid lines (cqs.|6]and[7]l. The uppermost 
panel shows the gas pressure evaluated as P = nT = 1.91neTe, 
where rie and Te are the best-fitting electron density and temper- 
ature. Outside the central 0'.4 the best fit parameters for the two 
bands agree well, suggesting that the impact of multi-temperature 
plasma on the rie and Tf, is not strong. 



Figure 6. 'Biases' in a two-temperature plasma versus volume 
fraction of the hotter component. Fake Chandra spectra repre- 
senting a mixture of fcTi = 1 and kT2 = 2 keV plasmas in pres- 
sure equilibrium were fitted with a one-temperature model in the 
0.6—9 and 2.5—9 keV bands (red and blue points respectively). The 
five panels from top to bottom show i) the best-fitting gas tem- 
perature, ii) the ratio of best-fitting gas density to the true mean 



density, iii) gas metallicity, iv) gas pressure Pfu 



T-fit-Lfit 



, v) 



the ratio of the measured value Pjit/nfit = Tfu to the true value 
P/rimean = Tmean] here Timean is the volume averaged plasma 
density (eq. [SJl . The last quantity characterizes the bias in mea- 
suring the gravitating mass or the potential from the hydrostatic 
equilibrium equation. 



the Galactic value Nh = 2.54 x 10^° cm"^ (Dickey & Lock- 
man 1990) and the redshift was fixed at the optically de- 
termined redshift of M87: z = 0.00436 (e.g., Smith et al. 
2000). The gas temperature, heavy metal abundance and 
normalization were free parameters of the model. A sam- 
ple spectrum in the 7'. 4 — 8'. 8 annulus and the deprojected 
spectrum in the spherical shell of the same size are shown 
in Fig. m together with the best fitting single-temperature 
APEC model. The best fitting temperatures are 2.57 and 
2.4 keV for projected and deprojected spectra respectively. 
The deprojected spectra were generated in narrow (14.6 eV) 
energy bins and used in XSPEC without any further group- 
ing. To avoid a bias towards lower temperatures, which is 
present when the number of counts per spectral bin at high 
energies becomes less than ~ 10 — 20, the statistical errors 
were evaluated from the smoothed spectra, while the fitting 
is applied to raw (unsmoothed) spectra (see Churazov et al. 
1996 for details). This approach suppresses the bias associ- 
ated with the low number of counts per spectral bin provided 
that total number of counts in the spectrum is large. We did 
verify earher for ASCA (Churazov et al. 1996) and XMM- 
Newton data and now for Chandra observations of M87 that 
this approach yields fully consistent results to an often used 
technique of grouping spectral channels in bins with at least 



10-20 counts (assuming that the correct spectral model is 
used) . Note that we also verified that allowing the Fe abun- 
dance to vary separately from the a-elements altered the fits 
by not more than the statistical errors. 

The electron density was derived from the normaliza- 
tion of the spectra, fixing the proton to electron ratio to 0.83 
and taking into account the distance to M87. The resulting 
spectral parameters are plotted in Fig. [5] as a function of 
distance from the center of M87. Red and blue points cor- 
respond to the best-fitting model parameters in the 0.6-9 
and 2.5-9 keV bands respectively. This choice of two en- 
ergy bands for spectral fitting is motivated on the following 
grounds: the broader band maximizes the statistical signifi- 
cance of the results, while the harder band provides a verifi- 
cation of the magnitude of biases which might arise when a 
one-temperature model is applied to a spectrum with mul- 
tiple temperature components (e.g., Buote 2000). 

To assess the magnitude of possible biases when fitting 
the multi-temperature plasma in the Chandra band with 
a one-temperature model, we generated a sequence of fake 
spectra representing a mixture of kTi — 1 and kT2 = 2 keV 
plasmas in pressure equilibrium, with a varying volume frac- 
tion /2 of the harder component. These two temperatures 
are characteristic for multi-temperature plasma in the inner 
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parts of M87 and NGC1399 (Buote 2002). The abundance 
of heavy elements was set to solar for both components, and 
Nh and redshift were fixed to the M87 values. The on-axis 
ACIS-I response was used. We then applied to the simu- 
lated spectra the same one-temperature model as we did 
for the real M87 data. The resulting spectral parameters 
are shown in Fig. [B] temperature, density, metallicity, pres- 
sure Pjit = UfitTfit, and the ratio of the measured value 



^\ ^ J ' ^*r 



Pfit/nfit = Tfit to the true value P/rir, 



-^ m 






ni X (1 -/2) 4-712 X /2 



(5) 



is the volume averaged plasma density and ni and 712 are 
the plasma densities of the cooler and hotter phases respec- 
tively (evaluated as 1.91 times the respective electron den- 
sities ris). The quantity Tfn/Tmean — 1 characterizes the 
bias in measuring gravitating mass or potential from the 
hydrostatic equilibrium equation. The red and blue points 
correspond to the best-fitting parameters derived from the 
0.6-9 and 2.5-9 keV bands respectively. For the 0.6-9 keV 
band, even a very small volume fraction of the cooler com- 
ponent causes the best-fitting temperature and metallicity 
to drop sharply, while the density rises (red points in Fig. 
[6|. For the 2.5-9 keV band, the density is, on the contrary, 
very close to the mean value. This conclusion is of course 
sensitive to the particular choice of Ti and T2 and to the 
effective area of the instrument. In general, one expects the 
best-fitting temperature in the broad band to be more bi- 
ased towards lower values compared to the hard band. It 
is possible therefore to use the differences in the best-fitting 
values obtained in two energy bands as an indicator of strong 
departures from a single-temperature plasma. From Fig. [5] 
it is clear that the best-fitting values in the 0.6-9 and 2.5-9 
keV bands differ strongly inside the inner 0'.4, while outside 
this radius, the results agree well. Of course even a reason- 
able agreement between the spectral results in two bands 
does not guarantee that the spectrum is fully described by 
a one-temperature model. This agreement nevertheless sug- 
gests that outside 0'.4 our spectral results are not strongly 
biased. On the other hand, the region inside 0'.4 is charac- 
terized by very strong departures from spherical symmetry 
and by the presence of cool and hot gas structures. Given the 
magnitude of the expected biases the results for this region 
should be taken with caution. 

From Chandra and XMM-Newton observations we 
know that the gas temperature in the X-ray bright features 
co-spatial with the extended radio bright 'arms' is lower than 
that of the typical intra-cluster medium (ICM) temperature 
at the same distance from the center (e.g., Belsole et al. 2001; 
Matsushita et al. 2002; Forman et al. 2005). These regions 
are clearly seen in Fig. [T] where the 0.6-2 keV image of M87 
is divided by an axisymmetric /3-model, as extended struc- 
tures to the east and the south-west of the nucleus. To test 
whether these structures dramatically affect our final result 
(e.g., the mass/potential profile, discussed below) we also 
fit the data excluding those regions containing the brightest 
parts of the arms. We use the spectral fits to this 'trimmed' 
data-set in the next section to verify the robustness of our 
results. 

Analytical approximations to the best-fitting values of 
ne{r) and kT{r) over the range of radii from ^ 0.5' to ~ 10' 





Figure 7. Central 10' X 10' part of the 0.6-2 keV image of 
M87, divided by the axisymmetric model with surface brightness 
I{R) oc (1 + (R/rc) 2)1/2-3/3 (j^gj.g ^ jg ^Yle distance from the 
M87 nucleus, /3 = 0.37 and re = 0.2') to show the structure in 
the surface brightness more clearly. 



l+(- 



/3 



are given below: 

He = 0.22 

where /? = 0.33 and r^ = 0.2' (or 0.93 kpc) 



kT= 1.55 X 



1 



\2.2'J 



keV. 



(6) 



(7) 



These approximations are shown in Fig. [S]by the black 
solid lines. 



5 POTENTIAL 

We first derive the gravitational potential for M87, including 
a detailed discussion of the method which we then apply to 
the data for NGC 1399. 



5.1 M87 

We assume that the gaseous atmosphere of M87 is in hydro- 
static equilibrium: 

I dP _ dip 
p dr dr 

where p is the gas density, P — P{r) is the pressure and (p 
is the gravitational potential. Thus 

1 dP 
p dr 

If the pressure is solely due to the thermal gas pressure, then 
P — nkT and p = prnpU, where p is the mean atomic weight 



V 



-dr 



(8) 



(9) 
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of the gas particles, nip is the proton mass, and n = n{r) is 
the total particle density. Then (assuming that fj, does not 
vary with radius) the potential tp is: 



Table 1. Data used to derive the potential. 



^ = 



/imp 



„dlnn , „ 
T^ — dr + T 
dr 



+ c 



(10) 



where C is an arbitrary constant. Below we will use the 
value of the potential relative to a given reference radius 
r^e/ = 5'. 5 (25.6 kpc), i.e., (y3(5'.5) = 0. Note that eq. (fllT)) 
depends only on the measured quantities n{r) and T{r). 
Furthermore, eq. (|10|l does not depend explicitly on radius 
r (only through n(r) and T{r)) or on the absolute normal- 
ization of the gas density. This in particular means that the 
assumed distance to M87 does not affect the calculated po- 
tential as long as it is expressed as a function of angular dis- 
tance from the nucleus. Thus <p can be easily evaluated using 
the densities and temperatures derived from the deprojected 
X-ray spectra. Unlike the expression for gravitating mass, 
which explicitly depends on the spatial derivatives of the 
observed quantities, the change of the potential A</5 between 
two radii ri and r2 depends primarily on the gas densities 
at two radii and on the mean gas temperature. For exam- 
ple, in the case of isothermal gas, A</5 oc Tin [n(r2)/n(ri)) 
which is of course a direct consequence of the Boltzmann 
distribution. This makes the calculation of Atp very simple 
and robust. 

Below we will denote the potential derived from X-ray 
data (by means of eq. I10|l as ipx to distinguish it from the 
true potential tp or from the potential tpopt derived from 
optical data. 

Note that the derived spectral parameters in each bin 
are not independent random variables since they come from 
the deprojection analysis, e.g., the emissivities in the spher- 
ical shells are various linear combinations of the same set 
of original independent data points (fluxes in the concentric 
annuli around the source). For a given spherical shell the er- 
rors are evaluated correctly since the spectrum is calculated 
as an explicit linear combination of independent variables. 
However, it is expected that the deviations from the correct 
values in nearby radial shells are correlated/anti-correlated. 
Thus, if we repeat the same observation many times, the 
measured values of the emissivity will lie within the esti- 
mated uncertainties in 67% of cases (apart from possible 
systematic errors) , but the values in nearby radial bins may 
change in correlated/anti-correlated ways. 

Similarly, the values of tpx obtained at different radii 
are not statistically independent, both because they are de- 
rived from the deprojected spectral parameters and because 
the value of tpx at each radius depends on the same set 
of spectral parameters in various radial bins according to 
equation (fTO|). Rigorous error calculations would require er- 
ror propagation from the original (independent) spectra ob- 
tained in the set of annuli around the source; these would be 
difficult to do and the results would be difficult to present 
simply. However, in practice we believe this more rigor- 
ous analysis may be replaced with a straightforward er- 
ror calculation. Indeed, for the case of an isothermal gas 



V3x(r) 



^ITip 



T In n(r) + C, i.e., ipx{r) depends only on the 



local (at radius r) gas density. Therefore, the error distri- 
bution for ifixir), at any given radius r, can be correctly 
estimated, as long as the errors in n{r) axe correct (see previ- 
ous paragraph) . The same is approximately true for a slowly 



Regions 



Energy band Color in Fig. [8] 



0-360° annuli 0.6-9 keV 

0-360° annuli, excluding 'arms' 0.6—9 keV 

0-360° annuli, excluding 'arms' 2.5—9 keV 



blue 
red 

green 



varying temperature profile. Below, we assign to each ipx (r) 
the error which follows from equation (|10l assuming that the 
errors in evaluating the spectral parameters at each radius 
are independent. 

In Fig. [8] we show three versions of ipx derived from the 
Chandra data for M87. Each of the three shaded 'curves' 
corresponds to different subsets of data and/or different en- 
ergy bands used for spectral fitting, as indicated in Table [l] 
The regions that exclude 'arms' are described more fully at 
the end of §4.2. The vertical width of each of the color curves 
in Fig. [8] corresponds to the 67% statistical uncertainty. 

All three methods lead to estimates of fxir) that are 
consistent with each other within statistical uncertainties. 
The potential derived from the data in the 2.5-9 keV band 
has the largest statistical errors, but as discussed above, 
the spectral parameters may be less biased in this band, 
compared to the broader 0.6-9 keV band. We note here that 
departures from spherical symmetry are obvious from Fig. [T] 
and we cannot expect the spectral parameters (or potential 
profiles) to be perfectly consistent with each other. 

For comparison to the X-ray derived gravitational po- 
tentials, we show in Fig. [S] with the black lines, two esti- 
mates of the gravitational potential that are derived from 
the optical data. The dashed line shows the potential de- 
rived by Wu & Tremaine (2006, their §4.3.2; WT06 be- 
low) using a sample of 161 globular clusters in M87. The 
mass profile in this model is a power law function of radius 
M{r) = 2.3 X 10" MQ(r/kpc)^=*®. This mass distribution 
is most accurately determined near r = 30 kpc, but Wu & 
Tremaine estimate that the relative error in M{r) is less 
than 40% for 17kpc < r < 90 kpc (or 3.7' < r < 19'). 
The error in the potential is substantially smaller, since 
it is determined by the integral of the mass distribution, 
(f^opt = J^ GM{r)dr/r^. The potential in the inner few 
arcminutes is only weakly constrained since this model does 
not use any information on stellar kinematics. The thick 
solid (black) line in Fig.|8]is the best-fitting model to the star 
and globular-cluster kinematics in M87 from Romanowsky 
& Kochanek (2001, see their §4.2; hereafter, RKOl). This 
model (NFW2 in their Table 2) has a mass profile that is 
the sum of an NFW profile (Navarro, Frenk & White, 1996) 
from a dark halo and a stellar component, derived from the 
M87 optical surface brightness under the assumption of a 
constant stellar mass-to-light ratio. In terms of the likeli- 
hood, this is the most probable model among those consid- 
ered by RKOl. This model agrees reasonably well with the 
X-ray derived potential apart from a few clear wiggles in 
ifix{r), which we discuss in Section [6.51 Two other models - 
NFWl and NFW3 (see Table 2 of RKOl) are shown in Fig. [8] 
with the thin solid lines. In terms of likelihood/significance 
these models are ~ 0.8 and 0.9 a "dispreferred" compared to 
the best fitting NFW2 model. The difference between these 
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Figure 8. Derived potential ipxix) using the 0.6-9 keV band 
(blue), the same excluding 'cool arms' (red), the same in the 2.5- 
9 keV band (green) — see Table[T] Black curves - potentials foptir) 
derived from the optical data: Wu & Tremaine (2006) - dashed 
line; Romanowsky & Kochanek (2001), NFW2 model - thick solid 
line. The comparison of the thick black and colored curves sug- 
gests that If opt (r) and ipxi''') broadly agree over the range of 
radii shown in the figure. Thin solid lines are the models NFWl 
and NFW3 from RKOl which in terms of likelihood/significance 
are ~ 0.8 and 0.9 a "dispreferred" compared to the best fitting 
NFW2 model. The difference between these curves approximately 
characterizes the uncertainties in ipopt (r) based on the analysis of 
RKOl. All curves are normalized to zero at 5'. 5 radius. 



curves can be used to approximately characterize the uncer- 
tainties in ipopt{r) based on the parametric model of RKOl. 
We use below the power law model of WT06 and the 
NFW2 model of RKOl as two versions of ipopt{r) and assume 
that the difference between these two models is comparable 
to the deviations of either from these models from true po- 
tential of the galaxy. 

5.2 NGC 1399 

The X-ray potential ifix{r) derived from Chandra data for 
NGC 1399 is shown in Fig. [9] with the blue curve. The po- 
tential is set to zero at l'.5. In X-rays the core of NGC 1399 
appears less disturbed than M87 and the resulting ifix (r) is 
much smoother than for M87. For comparison ipoptir) based 
on the model of Kronawitter et al. (2000) is shown with the 
black solid line, and their 2a confidence bands as the thin 
black lines. As in the case of M87, there is good agreement of 
ipx(r) and ipoptir) over the region where X-ray and optical 
data are available. 

The optical potential (jiopt for NGC 1399 was derived 
from absorption-line kinematic data, as described in Saglia 
et al. (2000). These authors constructed spherical models for 
the galaxy in a sequence of gravitational potentials, com- 



Figure 9. Derived potential ipxij) for NGC 1399 using the 0.6- 
9 kcV band for 0-360° annuli (blue curve). The thick black solid 
line shows the potential ipopt{r) derived from optical data of Kro- 
nawitter et al. (2000), while thin black lines show ±2(7 limits. Out- 
side the range of radii covered by the optical data (r < 1.62'), the 
extrapolation of the potential is shown by dashed lines. The green 
curve is the 'stars-only' model from Kronawitter et al. (2000). The 
magenta line shows the optical potential ipopt{i") multiplied by 0.9. 
All curves are normalized to zero at l'.5. The good agreement of 
fopti'i') and fxi''') suggests that the fraction of non-thermal pres- 
sure is low. 



posed of the potential of the luminous stars and various op- 
tional quasi-isothermal dark halo potentials. The potential 
of the stars was based on the deprojected luminosity distri- 
bution and the assumption of constant mass-to-light ratio. 
For each of these potentials, the non-parametric distribution 
function f{E,L^), that gave the best fit to the kinematic 
data (velocity dispersion and symmetric 4*'' moment /14 of 
the LOSVD), was found; here E and L denote energy and 
total angular momentum. The Ax^ values of the various 
models then determined the boundaries of the confidence 
region for the halo potential parameters. Similar models, 
based on a larger basis of distribution functions and poten- 
tials, were constructed by Kronawitter et al. (2000). The 
Kronawitter et al. models are plotted in Fig. |9] and span 
the 95% confidence range within their expanded set of po- 
tentials. The true confidence range based on non-parametric 
potentials is likely to be somewhat larger; however, it is not 
straightforward to determine. An independent investigation 
by Graham et al. (1998) yields similar results. 



6 EFFECTS OF MAGNETIC FIELDS, COSMIC 
RAYS, AND MICRO-TURBULENCE 

The good agreement between the gravitational potentials 
derived independently from the optical and X-ray data, 
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found in the previous section, suggests that hydrostatic equi- 
hbrium is satisfied reasonably well in M87 and NGC1399. 
This conclusion contrasts with that of Diehl & Statler (2007) 
who have argued recently that hydrostatic equilibrium of the 
X-ray gas in elliptical galaxies is "the exception rather than 
the rule" and "X-ray-derived radial mass profiles may be in 
error by factors of order unity". In particular, NGC1399 is 
one of the objects from the Diehl & Statler (2007) sample 
that is identified as being far from hydrostatic equilibrium 
based on differences in the isophotal shapes of the stars and 
the X-ray gas. We believe that the close agreement of the 
potential profiles 4>x{r) and (jioptir) obtained in the previ- 
ous section would be an improbable coincidence if hydro- 
static equilibrium were not approximately valid. While nei- 
ther NGC1399 nor M87 are in perfect equilibrium, as in- 
dicated by visibly disturbed X-ray images (see e.g. Buote 
& Tsai, 1995) , the estimate of the magnitude of the depar- 
ture from hydrostatic equilibrium by Diehl & Statler (2007) 
is too pessimistic. We will further test this statement in a 
larger sample of elliptical galaxies in future work. 

We now consider how magnetic fields, cosmic rays and 
micro-turbulence would affect the potential '^x{r) derived 
from the X-ray data. We keep the hydrostatic equilibrium 
assumption (eq. [5| , but we now allow the pressure or density 
derived from X-ray analysis to differ from the true P and p 
entering eq. ((9]). 

6.1 Uniform medium 

Let us first keep the assumption that the X-ray emitting gas 
is uniform within each spherical shell, but assume that this 
gas provides only part of the total pressure. The rest of the 
pressure is provided by a component that is 'invisible' in X- 
rays (e.g., cosmic rays), which is also uniformly distributed 
through the volume of each shell. Here and below we assume 
(unless stated otherwise) that the pressure is isotropic and 
the atmosphere is in hydrostatic equilibrium, with the ther- 
mal gas 'mechanically' coupled to other components. We 
define the fraction of the pressure due to thermal gas as 
fg = fg{r). Then the pressure derived from X-ray data will 
be nkT = fgP, where P is the true total pressure, n and T 
are the density and temperature measured from X-ray data. 
If we further assume that only thermal gas contributes to the 
mass density, then p = pnipn will be correctly determined 
from the X-ray analysis. The temperature is also assumed 
to be determined correctly. Substituting fgP instead of P 
into eq. (|9]), we can easily express <^x(r) through the true 
potential ip{r) and other gas parameters: 



<fix 



J iJ-mp J 



dr 



(11) 



If fg does not depend on radius then obviously ipx = fgf- 

Shown in Fig. [10] is ipx from eq. (|ll|l for fg = const — 
0.5 and 0.85 (magenta solid and dotted lines respectively) 
and for 



fg = [l + {r/l'f]/[2 + {r/l')' 



(12) 



(red dashed line). The latter case illustrates the situation 
when gas pressure dominates at large radii {fg = 1), but 
declines to 50% of the total pressure {fg = 0.5) in the core. 
In these calculations we used ip based on the NFW2 model of 
RKOl and the analytic approximation of temperature given 



by eq. ([7]). The curves for fg = l and fg = 0.85 fit the X-ray 
data reasonably well, while the other models are excluded. 
For NGC 1399, the ipx and (popt closely follow each other. 
Values of fg smaller than 0.9 (dotted magenta line in Fig. [9} 
would lead to a substantial disagreement between (fx and 
ifiopt over the range of radii from few arcseconds to ~ 2'. 



6.2 Micro-turbulence 

If small scale and isotropic turbulent motions are present 
in the medium, then the impact on ipx can be evaluated 
similar to the case considered in the previous section. The 
quantity 1 — fg characterizes the contribution of turbulent 
motions to the pressure (or energy density) of the medium. 



6.3 Bubbles of relativistic plasma 

We now assume that the medium is not uniform. Let fg 
be the volume fraction filled with thermal gas. The remain- 
ing volume fraction is completely devoid of thermal gas and 
is occupied by bubbles of cosmic rays and magnetic fields. 
Both phases are assumed to be in pressure equilibrium and 
coupled together. If, in performing the X-ray deprojection, 
one still assumes that the gas is uniform (as we did above), 
the derived gas density will be lower than the true gas den- 
sity in the gas patches and higher than the mean volume 
density; poba = Pgaas/Yg = Pm^anl\fTg- The pressure de- 
rived from the X-ray analysis will be modified accordingly. 
For this case, the potential is: 



<y3A- 



/grfV'- 



1 k 



2 /iTTlp 



dr 



(13) 



The X-ray potential <px derived from eq. (|13p for fg given 
by eq. (|12p is shown in Fig. [10] by the green line. 

We note here that a similar problem (the impact of 
bubbles on the observed gas distribution) has been recently 
considered in Nusser & Silk (2007). 



6.4 Summary of the effects of non-thermal 
pressure components on 'px{r) 

As discussed above we consider the following forms of non- 
thermal pressure support as the most relevant; (i) cosmic 
rays and magnetic fields uniformly mixed with the ther- 
mal gas (Section 16. If) , (ii) cosmic rays and magnetic fields 
forming bubbles that are free of thermal gas (Section I6.3|l 
and (iii) micro-turbulence in the thermal gas (Section 16. 2|) . 
For all these CEises non-thermal pressure manifests itself as 
a coefficient in the relation dpx ~ fgdp if fg is indepen- 
dent of radius. For fg varying with radius the dependence 
on fg{r) is slightly more complicated (see egs. 1111 and I13|) . 
with an additional term in (^x ('"2) — (l>x{ri) roughly equal 
kT , fg{r2) 



to 



■In^ 



,-, where fg{r2) and fg{r2) are the ther- 
prrip fg{ri) 

mal gas pressure fractions at two radii. Therefore, the com- 
parison of the change of ipx to the true potential (if it is 
known) over the broad range of radii from ri to r2 provides 
a simple and convenient way of estimating fg — the contribu- 
tion of the non-thermal components to the total pressure — 
independent of the "form" of the non-thermal pressure. This 
applies also to combinations of micro-turbulence, cosmic 
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Figure 10. Left: Various model dependences of the thermal gas fraction fg on radius discussed in Section[6] Right: Expected potential 
ifx that would be derived from X-ray observations for each of the fg models shown in the left plot. The color coding is the same in both 
plots. We assumed that the true potential is the NFW2 model of RKOl. Such a potential would be derived from perfect X-ray data if 
fg = const = 1 (thick black solid line). Other fg, ipx pairs shown on both plots are: solid magenta line represents fg = const = 0.5; 
dotted magenta line represents fg = const = 0.85; red dashed line represents fg given by eq. I I12II and fx from eq. ((TTJ; green solid line 
represents fg given by eq. II12I I and ipx from eq. 11131 . Note, that the gas pressure fraction curves for the two potentials in green and red 
(dashed) are identical and they overlap in the left panel. For comparison, the blue curve shows tfx derived from the Chandra observations 
of M87 (same curve as shown in Fig. [Sj. The comparison of black, magenta, red and green lines in the right figure demonstrates the 
impact of the non-thermal pressure on ipx- The measured ipx (blue line) agrees with the optical profile (black line) better than any of 
the model curves, suggesting that the fraction of nonthermal pressure is on average less than is assumed in any of these models. 



rays uniformly mixed with the thermal gas, and cosmic rays 
contained in bubbles. 



6.5 Sound •waves and shocks 

The presence of sound and shock waves formally invalidates 
our assumption of hydrostatic equilibrium. At any given lo- 
cation the potential derived from X-ray data can be either 
underestimated or overestimated, depending on the sign of 
the pressure gradient in the wave. Forman et al. (2005, 2007) 
identified in the X-ray data a number of quasi-spherical fea- 
tures around M87 which are plausibly caused by shock waves 
produced by an unsteady outflow of relativistic plasma from 
the black hole at the center of M87. These shock candidates 
are located exactly where we see wiggles in the potential 
(Fig, nop . This is not surprising since, in a spherical shock 
propagating from the center of the galaxy, the pressure first 
increases sharply inwards at the front and then decreases in 
the rarefaction region behind the shock. The potential de- 
rived from the X-ray analysis (see eq. llOf) will have a 'dip' 
at the position of the shock front and it will then recover 
from this dip over the rarefaction region. To illustrate such 
behavior we simulated a spherical shock in the atmosphere 
of M87 and calculated the potential ipx which would be de- 
rived from X-ray data corresponding to these simulations. 
Fig. llll shows the true potential, assumed in the simulations. 



with the black line, while the blue line shows ipx calculated 
for a Mach 1.2 shock located 2.7' from the nucleus (see For- 
man et al. 2007,2008 for details of the shock parameters and 
the simulations). 

For weak shocks/sound waves which are well localized 
in space (e.g., weak compression waves) px will be af- 
fected only in the region where the wave is present, while 
/S.ipx = ipx{r2) — ipx{ri) over a range of radii from ri to 
r2 which covers the entire perturbation will not be affected. 
Therefore, weak and localized perturbations do not affect 
the global change of px ■ In the discussions below that com- 
pare px and popt (Section 17. 1|) we focus on radii, ri to r2, 
away from the radii most strongly affected by the shock. 



6.6 Non-luminous gas 

In sections 16.11 16.2116.31 we have assumed that there is an 
'invisible' (in X-rays) component contributing to the pres- 
sure. As a result, the derived potential ipx was flatter than 
the true potential. We can also assume that there is an X- 
ray 'invisible' gas component which contributes to the mass 
density and has a small filling factor. A cool gas phase with 
temperature below 'observable' X-ray values but coupled to 
the X-ray emitting phase (e.g., by magnetic fields) would 
be the simplest example. The px derived from equation ((9)1 
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Figure 11. Potential {fx - blue line) derived using eq. IIIOI I from 
simulated density and temperature profiles. In these simulations a 
powerful explosion (with total energy a few lO^'^ erg) at the center 
of the galaxy drove a shock through the ICM. The data shown 
in the figure correspond to the time when the shock of Mach 
number ~ 1.2 is located ~ 2'. 7 from the nucleus, the observed 
radius of the strongest shock in M87. The true potential is shown 
with the black line. The deviations of (fix predicted by the shock 
model from the true underlying potential resemble the deviations 
between the observed tpx and (fiopt (e.g. blue curve in Fig llOI I. 
suggesting that most of the wiggles seen in the observed ipx arc 
due to the shock. 



will be simply: 
Vx = frndip, 



(14) 



where fm ~ p/px is the ratio of total mass density to the 
mass density derived from the X-ray analysis. In this case 
fm > 1 and the X-ray derived potential will be steeper than 
the true potential. Although we note this possibility here, 
for the remainder of the discussion below, we assume that 
this "non-luminous gas" component does not affect M87 and 
NGC1399. Otherwise it will be difficult to disentangle the 
contribution of such a component from that of non-thermal 
components which affect ifix in the opposite direction. In 
the absence of any "non-luminous gas" component, we can 
place constraints on the non-thermal pressure as we discuss 
below. 



7 DISCUSSION 

7.1 Evaluation of the non-thermal pressure 
contribution fg 

In Fig. [12] the potential derived from X-ray data is plot- 
ted against the potential derived from optical data for M87 
and NGC 1399 (red shaded areas) . If the assumptions used 



when deriving the potentials are correct, then one expects 
to find ipx = ipopt (shown with the dashed lines). If cos- 
mic rays, magnetic fields or micro-turbulence contribute to 
the pressure in the X-ray emitting gas, then the change of 
ifix between two radii will be smaller than the change in 
fopt, i.e. Aifix = aAipopt, where a < 1. In the simplest 
case of a constant fg, there is a linear relation between the 
two potentials tpx = afopt + b, where a = fg. As shown in 
sections 16.11 and 16.31 if fg depends on radius then the rela- 
tion between the potentials can be more complicated. If fg 
is not constant, but decreases towards small radii, i.e. the 
role of non-thermal pressure increases towards the center, 
then ifix will be even shallower than is prescribed by the 
relation tpx = fgfopt + b (see eq. 1111 and I13|l . Therefore, 
wo can use the relation ipx ~ aipopt + b, corresponding to 
the assumption of constant fg = a, to provide a conserva- 
tive upper limit on the contribution of non-thermal pressure 
1 — fg over the observed range of radii. This upper limit cor- 
responds to an averaged non-thermal pressure contribution 
over the observed range of radii. The averaging is done with 
effective "weights" corresponding to the values of the opti- 
cal potential at each radius. For example, for an isothermal 
(logarithmic) potential the weight of the radial range from 

, . . , , r -I- Ar 

r up to r -f- Ar is proportional to In . 

r 

A linear fit to the observations of the form (fix = a'Popt + 
b is shown in Fig.[l2]with the thick solid line. For NGC1399, 
a = 0.93. For M87 a = 0.89 if the WT06 model is used and 
a = 0.975 for the NFW2 model of RKOl. The very closeness 
of a to unity in the latter case (NWF2 model) is surely a 
coincidence given the obvious wiggles in ifixir) (Fig. I12p . 
For the two other (less probable) models of RKOl - NFWl 
and NFW2 the values of a are 0.96 and 1.21 respectively. 
By repeating the fit four times using independent quadrants 
of the X-ray data (see Section 7.2) we conclude that the 
statistical uncertainty in these slopes is < 0.04. Taken at face 
value, these results mean that the fractional non-thermal 
contribution to the pressure amounts to 7% for NGC1399 
and 2.5-11% (depending on the ipopt model) for M87. These 
values characterize the non-thermal contribution averaged 
over the range of radii from 0'.7 to 12' in M87 and from 
3" to 5' in NGC 1399. Given that systematic uncertainties 
are present in the data (e.g., due to the shock in M87) it 
is difficult to provide a precise estimate of the uncertainty 
associated with this value. If one adopts a more conservative 
approach of upper limits (rather than measurements) on the 
non-thermal pressure components, then these limits are ~ 
10% and ~ 20% for NGC 1399 and M87 respectively (see 
Fig. [121). 

Fig. |6] (top panel) shows that, for temperatures char- 
acteristic of M87 and NGC1399, fitting multi-phase plasma 
spectra with a single-temperature model will result in (px 
biased towards shallower values (smaller a). Therefore this 
"multi-phase bias" has the same sign as the non-thermal 
pressure and can only lead to an overestimation of the non- 
thermal pressure fraction; thus the upper limits given in the 
preceding paragraph remain valid. In addition, based on the 
comparison of spectral fitting in two energy bands (Fig. [S|, 
we do not expect this bias to be strong, at least outside the 
innermost 0'.4 in M87. 

Errors in the optical potential profile could also con- 
tribute to the observed difference in potentials. In both 
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galaxies, uncertainties arise from the assumption of spheri- 
cal symmetry, the restriction to a particular family of para- 
metric halo potentials in the dynamical modeling, and pos- 
sible systematic errors in deriving the line-of-sight velocity 
distribution from the observed spectra. These uncertainties 
are difficult to estimate precisely. Confidence limits within 
the parametric potentials used for NGC 1399 are given in 
Fig. [9] and Section [5]2l Thomas et al. (2007) have inves- 
tigated the uncertainties in the inferred circular velocity 
curves for a sample of Coma cluster ellipticals which they 
studied with axisymmetric dynamical models and in a wider 
range of potentials. In these galaxies, the central mass den- 
sity is dominated by the luminous matter, independent of 
the dark matter halo profile used; thus the circular veloci- 
ties inside ~ 0.5-Re are accurately determined, within ~ 5%. 
At two effective radii, the uncertainties in the inferred circu- 
lar velocities for their galaxies are ~ 15%. The measurement 
errors in the observed kinematics for NGC 1399 (Saglia et 
al. 2000) are smaller than in the Coma ellipticals. However, 
NGC1399 observations are along only one long slit. Hence, 
we would expect somewhat smaller uncertainties in NGC 
1399. Again the error in the potential is smaller, since the 
potential is obtained by the integral of the circular velocity 
curve. However, note that differences (f)opt — (jitrue, in the in- 
ferred optical potential profiles are likely to be more slowly 
varying functions of radius than the variations in (f>x — (jiopt 
seen in Fig. O 

The ifx — fopt difference characterizes the combined 
contribution of all non-thermal pressure components and the 
modelling uncertainties of the respective datasets. Under the 
assumption that modelling uncertainties do not dominate 
and that the contributions of individual components of the 
non-thermal pressure all have the same sign, the measured 
values can be converted to upper limits. For example, frac- 
tions of non-thermal pressure support of 10% in M87 and 
NGC 1399 translate into upper limits on the magnetic field 
of ~16 and ~10 fj,G respectively (evaluated at a distance 
from the galaxy center of 2' in M87 and 1' in NGC 1399). 

7.2 Non-spherical models 

We have assumed that the gravitational potential in both 
M87 and NGC 1399 is spherical. The eUipticity of the opti- 
cal isophotes of M87 is near zero at the center, rising to 0.4 
at 10', while the eUipticity of the X-ray isophotes outside 5' 
is 0.1-0.15, with the same position angle (Bohringer et al. 
1997) (see Fig[T} . The eUipticity of the optical isophotes of 
NGC 1399 is about 0.1 out to 10' (Dirsch et al. 2003). The 
X-ray image of NGC1399 looks more disturbed (see Fig. [2]) 
than the optical image and has an eUipticity of 0.34 at the 
effective radius of the galaxy (Diehl & Statler 2007). The 
shapes of the X-ray emitting gas and the stellar distribu- 
tion need not be the same in the common potential, since 
the velocity-dispersion tensor of the stars is not necessarily 
isotropic. 

Thus both galaxies exhibit a modest level of non- 
axisymmetry in their X-ray and optical isophotes. We can 
assess the impact of these asymmetries on the X-ray de- 
rived potential by using the data in individual (indepen- 
dent) wedges and deriving the corresponding X-ray poten- 
tials as shown in Fig. [I3l Within each 90° wedge, cen- 
tered at the optical center of the galaxy, we assume that 



the gas properties depend only on the distance from the 
center and we repeat our previous analysis up to fitting 
the (fix = atfiopt + b relation. The innermost and outer- 
most bins, which have large statistical uncertainties, were 
excluded from the fit (as shown by thin vertical lines in 
Fig.[T3}. The values aM87,NWF2 = {1.04, 0.98, 1.03, 0.95} and 
aNGCi399 = {0.95, 0.93, 0.92, 0.95} were found for the NW, 
NE, SE and SW wedges in the two galaxies. If WT06 model 
is used for M87 then aM87,WT06 = {0.93,0.86,0.89,0.87}. 
The root-mean-square deviations of a from the mean values 
are 3-4% for M87 (depending on the optical profile used) 
and 1.5% for NGC1399. These modest deviations character- 
ize the total uncertainty associated with the apparent non- 
axisymmetry in the X-ray isophotes of these two galaxies. 
Note that the observed scatter in slope between the four 
quadrants also represents an upper limit to the statistical 
uncertainties in the slope, although for these datasets we 
believe that the scatter is mostly due to non-axisymmetry 
rather than statistics. 

An additional uncertainty, probably of a few per cent, 
is related to the freedom of choosing the range of radii used 
for fitting the (px = afopt + b relation. This is more serious 
for M87 where the presence of a shock clearly makes the 
results more sensitive to the choice of the radial bins. 

Another concern is that the actual deviations from 
spherical symmetry may be larger than the observations in- 
dicate, if the galaxy is an approximately axisymmetric sys- 
tem that is viewed pole-on. The effects of a non-spherical 
potential on X-ray estimates of mass or potential have been 
modelled by Piffaretti, Jetzer & Schindler (2003) for an el- 
lipsoidal /3-model (see also Gavazzi 2005 for the analysis of 
axisymmetric NFW profiles) . Although Piffaretti et al. were 
concerned with clusters of galaxies, their results apply to 
galaxies as well. They found that if (i) the galaxy potential 
is axisymmetric and viewed pole-on and (ii) its mass is es- 
timated from the temperature and distribution of hot gas 
with the assumption that the potential is spherical; then, 
the derived mass within a spherical radius is typically in 
error by only a few per cent, if the axis ratio of the poten- 
tial is varied within plausible limits (typically ±40%). We 
note that for the model considered by Piffaretti et al. (2003) 
(ellispoidal /5-model distribution of isothermal gas) in the 
limiting case of the radius much larger than the /3-model 
core radii the deprojected potential/mass profiles are exact 
for any orientation of the galaxy/cluster with respect to the 
line of sight (see Appendix[B]| . While the above assumptions 
(e.g. isothermality) are not in general met, this result sug- 
gests that the effect of triaxiality can be small. We provide a 
more quantitative analysis of the triaxiality on real systems 
in future work. 

The effects of a non-spherical potential on optical esti- 
mates of mass are more difficult to model. The simplest ap- 
proach is based on the tensor virial theorem (e.g., Binney & 
Tremaine 2008). If (i) an axisymmetric one-component self- 
gravitating stellar system is viewed pole-on; (ii) the density 
of stars in the system is constant on spheroids of eccentric- 
ity e; and (iii) the mass of the system is estimated using the 
virial theorem for spherical systems; then the derived mass 
will exceed the actual mass by a factor /, where 
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Figure 12. Potential derived from X-ray data plotted versus tlie potential derived from optical data for M87 and NGC 1399 (red shaded 
areas). The thick solid line is a formal linear fit tfx = dfopt + b. In M87 the value of a is 0.975 for NWF2 model of RKOl and 0.89 
for WT06 model. In NGC 1399 a = 0.93 for Kronawitter et al. 2000 model. The thick dashed line shows the dependence i^x = Vopti 
while thin dotted lines correspond to a=0.8 and 1.2. For both objects, values of a close to unity imply that the fraction of non-thermal 
pressure is ^0.1. 
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of the potential in different ways. The close agreement be- 
/, r-)tween the X-ray and optical potentials, seen in Figs. |5]and 



would be an unlikely coincidence if any asphericity were 
significantly affecting the derived potentials. 



where e = ^1 — W- jd^ and hja is the axis ratio. For a plau- 
sible range of axis ratio, say 0.5 < hja < 1, / varies only 
between 0.79 (most oblate) and 1.04 (most prolate), suggest- 
ing that the error in the optical potential due to asphericity 
is likely to be small. 

It should also be stressed that the probability that ei- 
ther of the studied galaxies is strongly oblate or prolate and 
viewed pole-on is rather low. The distribution of ellipticities 
in a sample of brightest cluster ellipticals similar to M87 and 
NGC1399 was studied by Porter et al. (1991). They mod- 
elled the distribution of axis ratios q as a Gaussian with 
mean \i, and dispersion (Je, and found /i = 0.66, a^ — 0.09 
if the galaxies were assumed to be oblate, and /i = 0.69, 
(Te = 0.09 if the galaxies were prolate (at a radius of 15 kpc, 
corresponding to 3'. 2 in M87 and 2'. 6 in NGC 1399). The 
apparent axis ratios in M87 and NGC 1399 at this radius are 
0.2 and 0.1 respectively. The conditional probability distri- 
bution of the intrinsic axis ratio, given the apparent axis 
ratio, can be computed from these data. The probability 
that the intrinsic axis ratio in the stellar light distribution 
is as small as 0.5 is only 2% if the galaxies are oblate and 
6% if they are prolate. The potential due to the stars is, of 
course, more nearly spherical than the distribution of the 
stars themselves. 

Finally, we point out that substantial errors due to as- 
phericity are likely to affect the optical and X-ray estimates 



7.3 Is the small non-thermal component 

consistent with AGN-heated cluster cool 
cores? 

The central region of M87 is a classic example of a cooling 
flow - a region where gas radiative losses are high and an 
external source of energy is required if one wants to have a 
quasi-steady state without large net cooling of the gas, as 
indicated by the recent XMM-Newton RGS data (e.g., Pe- 
terson et al. 2003). A plausible source of energy for gas re- 
heating is the activity of the central black hole, mediated by 
the outflow of relativistic plasma. Actual gas heating might 
be due to dissipation at shock fronts (e.g., David et al. 2001) 
or in sound waves (if the medium is viscous, has apprecia- 
ble thermal conductivity or is strongly non-uniform, e.g., 
Fabian 2003, 2005; Ruszkowski et al. 2004a, 2004b; Heinz & 
Churazov 2005) or through an intermediate step of gener- 
ating gas turbulence, which eventually dissipates its energy 
(e.g., Churazov et al. 2001, 2002). We now test whether our 
results on the non-thermal forms of pressure are consistent 
with this picture. 

First, we directly see bubbles of relativistic plasma (as 
radio bright regions and as 'cavities' in X-rays). If the M87 
gas is in a quasi-equilibrium state with cooling losses, on av- 
erage compensated by heating, then the mechanical power of 
the outflow should be related to the gas cooling losses. As ar- 
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Figure 13. Potentials derived from X-ray data for M87 (left) and NGC1399 (right) in four 90° wedges (green, magenta, cyan and red 
for NW, NE, SE and SW wedges respectively), plotted against the optically derived potentials (RKOl for M87 and Kronawitter et al. 
2000 for NGC1399). The centers of the wedges coincide with the optical centers of the galaxies. The blue curves in each plot show the 
X-ray potentials derived from the full 0-360° data. Black solid lines show the linear least-square fits ipx = o-'^opt + b for each wedge and 
for the complete dataset. The innermost and outermost bins, which have large statistical uncertainties, were excluded from the fit (thin 
vertical lines). As in Fig. 1121 the thick dashed line shows the dependence <^x = Vopt, i.e. a = 1, while the thin dotted lines correspond 
to a=0.8 and 1.2. Clearly the slopes a derived for individual independent wedges agree well with each other, with root-mean-square 
deviations of 4.5% for M87 and 1.5% for NGC1399. The deviations characterize the overall uncertainties in a that are caused by the 
apparent asymmetries in the X-ray images and by statistical errors. 



gued in Churazov et al. (2001), Begelman (2001), and Nulsen 
et al. (2006), the efficiency of energy transfer from a rela- 
tivistic outflow to the gas is high, when the outflow subson- 
ically crosses regions with large pressure gradients. In this 
case, one can simply assume that the power of the outflow is 
equal to the cooling losses. The total energy of the relativis- 
tic plasma bubbles in the region can then be estimated from 
the balance of powers as E/tcroas ~ Sthermai/tcooi, where 
^thermal IS the gas thermal energy, tcooi is the gas cool- 
ing time and t cross is the time each bubble spends inside 
the region where it deposits a significant fraction of its en- 
thalpy. For example, consider a 10 kpc region around M87. 
For an electron density rie = 2 x 10~^ cm~'^ and a temper- 
ature Te = 1.7 keV, the cooling time tcooi ~ 7 x 10* yr. 
On the other hand, the rise velocity of large buoyantly 
driven bubbles is cCs where Cs is the sound speed and <; 
is a factor of order 0.5 (e.g., Churazov et al. 2001). Thus, 
tcross « 3 X 10^ yr for the 10 kpc region. Therefore, bub- 
bles of relativistic plasma should, on average, account for a 
fraction tcross/tcooi ~ 0.05 of the total energy (or volume). 
This is consistent within a factor of 2 with the results we 
obtained from Fig. 1121 Note that since the electron density 
in a typical cooling flow (e.g., in M87) scales approximately 
as We oc r~^ , then tcooi oc l/uc ex r and tcross oc r. Therefore 
the above estimate of tcross /tcooi is only weakly sensitive to 
the radius used to evaluate all quantities. 

We can now check whether our limits exclude turbu- 



lence as an intermediate reservoir of energy in the process: 
moving bubbles — » turbulence — » viscous dissipation and gas 
heating as suggested by Chandran (2004) and Rebusco et al. 
(2005, 2006). In this case, the energy dissipation rate of the 
turbulent motions should be approximately equal to the gas 
cooling losses: 
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where vturh and hurh are the characteristic velocity and spa- 
tial scale of energy containing eddies, n is the gas particle 
density and C is a dimensionless constant of the order of 
unity. According to a recent compilation by Dennis & Chan- 
dran (2005), C ~ 0.42. If the energy in the micro-turbulence 
is constrained to be less than a fraction ^ of the thermal en- 
ergy, then 



and then 



Vturb ^ 11 1 
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km s . This con- 



1 kpc^ 

dition is likely satisfied in the gas around M87 (if it is indeed 
turbulent), since the characteristic size of individual bubbles 
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is of order 1-3 kpc and their velocities are expected to be 
significantly higher, a few 10^ km s~^ (Churazov et al. 2001). 

7.4 The history of the ICM enrichment with 
cosmic rays 

Our constraints on the non-thermal pressure are consistent 
with the current paradigm of re-heating of gas in "cooling 
flows", which assumes a balance between radiative cooling 
and mechanical heating of the gas. On the other hand, our 
conclusion that the fraction of the energy density in the form 
of cosmic rays and magnetic fields is < 10-20% provides an 
important constraint on the history of the ICM. In particu- 
lar, these limits apply to the mixing of thermal plasma with 
relativistic plasmas, as well as the generation of cosmic rays 
by shocks in the bulk of the ICM. Generally, it is believed 
that the energy density of cosmic rays is dominated by pro- 
tons rather than electrons and the particles with the lowest 
7 provide the most important contribution to the energy 
density. The lifetime of a transrelativistic proton (Lorentz 
factor of 2) with respect to Coulomb energy losses in a clus- 
ter plasma with electron density of 10~^ cm^^ is roughly 5 
Gyr. Therefore, the energy density of relativistic protons is 
largely conserved, except for adiabatic losses, during the life- 
time of the cluster (transrelativistic electrons have of course 
much shorter lifetimes). This implies that our constraints on 
the energy density of cosmic rays in the cores of M87 and 
NGC 1399 today actually limit the presence of relativistic 
protons during the entire history of the same gas lump, pro- 
vided that relativistic particles are not able to diffuse out of 
the lump. The constraints on cosmic ray energy density in 
the past may be even tighter, if cosmic rays were originally 
deposited into a gas lump which had much higher entropy 
than is observed today and has lost the excess entropy by 
cooling. Indeed, if relativistic protons experience only adi- 
abatic losses, then their energy density changes with gas 
density as: 
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where pgas and £cr denote the thermal gas density and the 
energy density of cosmic rays observed now; 7cr = 4/3 is 
the adiabatic index of cosmic rays. All quantities measured 
at a different time are denoted with a prime. The thermal 
energy is by definition 
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Therefore, a constraint on the energy density of cosmic rays 
now Ecr/sgas ^ 0.1, for M87 and NGC 1399, can be trans- 
formed into constraints on e'cr/s'gas in the same gas lump. 
This is illustrated in Fig.[l4l Solid and open squares show the 
gas density and temperature observed in the cores of M87 
and NGC 1399 today. Along the solid lines (T oc p^"'"'^), the 
ratio e'cr/e'gas is constant. Therefore, if, in the past, the gas 
lump we see now in the core of M87 or NGC 1399 had a den- 
sity and a temperature along the thick solid line, then the 
constraint on the energy density of cosmic rays would be the 
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Figure 14. Constraints on the energy density of cosmic rays 
during the evolution of gas lumps presently seen in the cores of 
M87 and NGC 1399. The ratio c'cr/^'gas i^ constant along the solid 
lines. If the gas had in the past a density and temperature along 
the thick solid line, then at that time the constraint on energy 
density of cosmic rays is the same as now: e'^r/^'gas ^ O-I- Thin 
solid lines show the parameters for which present day constraints 
transform into e'cr/^'gaa limits of 0.01 and 1, respectively. For 
comparison we show a few typical states of the IGM/ICM, ranging 
from warm gas in sheets and filaments (triangle symbols and the 
dash-dotted line) in the large scale structure to shock-heated gas 
(3 and 8 keV) in galaxy clusters (dotted curves). The dashed line 
shows the locus of (rie , Te ) points taken from radially averaged 
profiles of the ICM in M87 at distances from 0'.5 to 5' from the 
galaxy center. 



same as now: e'cr/s'gas ^ 0.1. The two thin solid lines show 
the parameters for which present day constraints transform 
into Ecr/Sgas fimits of 0.01 and 1 respectively. For compari- 
son, we show a few typical states of the IGM/ICM, ranging 
from warm gas in sheets and filaments (triangle symbols 
and the dash-dotted line) in the Large Scale Structure to 
shock heated gas (3 and 8 keV) in galaxy clusters (dotted 
curves) . The cluster data are taken from simulations of Na- 
gai, Vikhlinin & Kravtsov (2007) for z = Q. The data for 
3 keV and 8 keV clusters were used and the radial profiles 
plotted in Fig. 114! correspond to distances from ~70 kpc to 
~2 Mpc from the cluster center. The main conclusion of 
this simple exercise is that if the gas in the core of M87 or 
NGC 1399 were in the past shock heated above the upper 
solid line, the energy density of cosmic rays in the gas, at 
that time, was less than 1-2% of the thermal energy density. 
We emphasize again that this conclusion is based on the as- 
sumption that cosmic rays dominating the energy density 
(presumably low energy protons) are not able to diffuse out 
of the gas lump and they suffer only adiabatic losses. Allow- 
ing diffusion would essentially imply that our upper limits 
are applicable only to the moment of observation. 
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Gamma-ray observations of clusters (Reimer et al. 2003) 
do limit the presence of cosmic rays, in particular protons, 
to less than ~20% of the gas thermal energy (e.g., Pfrommer 
& Ensshn 2004). It is interesting that Brunetti et al. (2007) 
have recently derived an upper limit on the contribution of 
relativistic protons to the gas pressure based on radio up- 
per limits for clusters without radio halos. While the actual 
limit on protons is sensitive to the assumed strength of the 
magnetic field, for many clusters the upper limit falls below 
1% (for > jiG fields) of the thermal pressure measured from 
the X-ray gas. 



7.5 Clusters of galaxies 

X-rays can be used to constrain mass/potential profiles of 
elliptical galaxies as long as the galaxy is gas rich and the 
contribution of individual stellar-mass systems to the X-ray 
emission (e.g., Revnivtsev et al. 2007) can be neglected. This 
implies that the X-ray analysis described in this paper can 
be applied to massive (> 10^^ Mq) systems. The most mas- 
sive ellipticals often dwell at the centers of groups/clusters 
(galaxies hke NGC1399 and M87). For these systems, X- 
ray observations can be used over a broad range of radii, 
from a fraction of the central galaxy effective radius to dis- 
tances almost comparable to the cluster virial radius. The 
optical data on stellar kinematics are mostly useful up to 
1-2 Re- This range of radii can be extended to a few Re 
with tracers such as globular clusters and planetary nebula 
(e.g., Samurovic & Danziger 2006, Wu & Tremaine 2006, 
Douglas et al., 2007), and to much larger radii using satel- 
lite galaxies or weak gravitational lensing. The comparison 
of lensing and X-ray masses made by Mahdavi et al. (2007) 
for a sample of 18 clusters yielded Mx/Ml = 1.03 ± 0.07 
and Mx/Ml = 0.78 ± 0.09 at radii corresponding to over- 
densities of 2500 and 500 respectively. These results, which 
characterize cluster properties at radii ~0.3-1.5 Mpc, are 
broadly consistent with our findings for the inner few tens 
of kpc. When combined with our work, these results sug- 
gest that the contribution of non-thermal components to 
pressure does not exceed ~10-20% over the entire cluster 
volume, as long as the system has not recently experienced 
a major merger. In numerical simulations of structure for- 
mation there is a substantial (in the range from a few per 
cent to tens of percent) contribution from residual stochas- 
tic gas motions to the pressure (e.g., Evrard 1990; Nagai, 
Vikhlinin & Kravtsov 2007), especially in the outer regions 
of clusters. If one accounts for this apparently unavoidable 
contribution, the remaining room for cosmic rays and mag- 
netic fields shrinks further, to ^10% over the whole cluster 
volume. Note however that the conditions in the cluster out- 
skirts and core are markedly different and it is unlikely that 
a universally applicable limit can be easily derived. 



8 CONCLUSIONS 

Using Chandra observations of M87 and NGC 1399 we de- 
rive the gravitational potential profiles of these two galaxies 
and compare them with the potentials derived from opti- 
cal data. Both the X-ray and optical methods - one based 
on the hydrostatic equilibrium of hot gas and the other on 
stellar dynamics - give consistent changes of the potential 



to within ^^10-20% over a broad range in radial distance. 
In M87, several obvious wiggles in the X-ray derived profile 
can plausibly be attributed to a spherical shock propagating 
through the ICM. 

These results suggest that the combined contribution 
of cosmic rays, magnetic fields and micro-turbulence < 10- 
20% of the gas thermal energy. These findings are consistent 
with the current paradigm of cool cluster cores, based on the 
assumption that AGN activity regulates the thermal state 
of the gas by injecting energy into the ICM. 

We also show that the limit of ~ 10-20% on the energy 
density in the form of relativistic protons applies not only 
to the current state of the gas, but essentially to the entire 
history of the ICM, provided that cosmic ray protons evolve 
adiabatically and their spatial difi'usion is suppressed. 

The comparison of tpx and ipopt is presented here for two 
objects only. Given the uncertainties in optical and X-ray 
modelling and peculiarities of the considered objects (e.g. 
notable shock in M87) it is desirable to extend the analysis 
to a larger sample. This project is under way and the results 
will be reported in subsequent publications. 
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APPENDIX A: A NOTE ON DEPROJECTION 
TECHNIQUES 

Mathematically the deprojection procedure described in 
Section 3] is identical (except for the l/rj factor) to the con- 
ventional onion-peeling algorithm (e.g., Fabian et al. 1980) 
if the number of spherical shells is equal to the number of ra- 
dial annuli and a unique and exact solution exists. However 
our formulation (i.e. eq.[3J works also for diff'erent numbers 
of shells/annuli. Typically Us < Ua, so the inversion prob- 
lem is overdetermined, but problems with tIs > Ua are also 
easily tractable with the additional condition that the solu- 
tion should have a minimal norm. Our method also provides 
an explicit expression for the statistical uncertainties in the 
solution. 

Compared to XSPEC's pro jet our deprojection 
method is much faster, given the large number of radial 
shells used here, since it avoids simultaneous fitting of many 
parameters. The expectation values are the same for our 
method, the onion-peel algorithm and XSPEC's pro jet if 
the spectral model is the correct one and the assumption 
of spherical symmetry holds. If, however, one of the above 
assumptions is incorrect, none of the approaches yields a cor- 
rect result. The proj ct results can be made more physical by 
imposing constraints on the parameters (e.g., non-negative 
fiuxes in every spherical shell), but at the expense of solving 
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a much more difficult problem of finding the true minimum 
in a multi-dimensional space. 

One particular application of XSPEC's the onion peel- 
ing pro jet that is sometimes used is fitting one spherical 
shell at a time and moving from outside in. This is faster 
since only one spectrum is fitted at each step. This proce- 
dure requires the spectral model to be specified in each outer 
shell to be able to fit the inner shells, since the contribution 
of these shells to the inner shells is evaluated using the spec- 
tral model (convolved with the instrument effective area and 
resolution). In contrast in our approach the observed I/77 
spectra are "subtracted" from the inner shells (after appro- 
priate scaling to account for geometry of shells) and the 
spectral analysis deferred to the final step of fitting depro- 
jected spectra. Hence, any incorrect model assumptions are 
not propagated into all the interior shells. Therefore, our 
deprojection analysis requires only assumptions about the 
geometry of the cluster. 

The bottom line is that all techniques should yield con- 
sistent results if underlying assumptions (e.g. spherical sym- 
metry and spectral model) are correct. Our depojection pro- 
cedure is especially simple and efficient if many radial bins 
(spherical shells) are used. 



APPENDIX B: X-RAY MASS AND POTENTIAL 
OF A NON-SPHERICAL CLUSTER. 

Here we prove the following result, mentioned in Section r7.2l 
Consider gas that is in hydrostatic equilibrium in a non- 
spherical gravitational potential. Assume that (i) the gas 
density is characterized by a power law in radius (a large 
radius limit of the ellipsoidal /3- model) , that is. 



p,as{r,e,4>)^f{e,4>)r- 



(Bl) 



where (r, 6, (j)) are the usual spherical coordinates and f{6, (f)) 
is an arbitrary positive function; and (ii) the gas is isother- 
mal. Then the gravitational mass within a sphere of given 
radius and the difference in gravitational potential between 
any two radii at the same {9, <f)) axe given correctly by an 
analysis that assumes that the gravitational potential and 
the gas density are spherically symmetric. 

To prove this, we write the hydrostatic equilibrium 
equation for an isothermal gas, 

Pgas{r,6,(t>) (X ex.p[-iimpip{r)/kT]. (B2) 

Equating this to HB1|I and taking the log, we have 



<pir,e,4>) 



kT 
p,mp 



[7 log r - log f{e, (t>)] . 



(B3) 



The corresponding density distribution is 



p{r,e,(k) 



47rG 
kT 



VV(r,e,0) 



[J-v'inm 



(B4) 



47rG/i77ip Lr^ 

Here we implicitly assume that f{0, (f>) corresponds to a 
physically sound density distribution, i.e., p{r,9,(f>) ^ in 
eq. HB4|) . For this distribution the gravitating mass within a 
sphere with radius r is: 



M{R) 



drp{r,e,4)) 



kT 
finipG 

-ykT 
l-tmpG 



7r - 



/ 

^|r|<fl 



drV^lnf{9,,p) 



(B5) 



since from the divergence theorem J. , „drV^ln/(6',0) = 
f. ,_ V In f{0,(t>) ■ ds — 0, where the second integration is 
done over a surface of a sphere with radius R. 

If we assume spherical symmetry when analyzing this 
cluster (with arbitrary orientation of the main axis relative 
to the line of sight), we will of course recover the same power 
law dependence of the density, px oc r~^. Then the change 
in the potential derived from X-ray data will be: 

Aipx = Vx{r2) - 'px{ri) 

^^ In pxin)/ px{r2) ^ ^^Hr2/ri). (B6) 



prrip 



pnip 



This is exactly the change of the true potential along any 
direction (eq. IB3|I . 

The derived mass within radius r will be: 



Mx{r) 



1 dPx r^ 1 kT d log PA- 



G px dr G r fiTUp dr 



-ykT 
l-impG 



r, (B7) 



where Px ~ kTpx / pnip. This is the correct value of the 
mass within the sphere of radius r (eq. |B5j. Thus, if the 
above conditions are met (isothermal gas with a power law 
density dependence), the X-ray analysis yields the correct 
spherically averaged mass and potential difference along any 
radius vector, independent of the orientation of the triaxial 
potential with respect to the line of sight. In other words 
- treating the triaxial system as spherically symmetric and 
performing a deprojection analysis as for a perfectly spheri- 
cally symmetric cluster yields the correct answer on the po- 
tential/mass. The same is true for the X-ray analysis applied 
to any wedge, no matter how narrow it is. 
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